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ABSTRACT 

We present new weak lensing observations of 1E0657--558 (z — 0.296), a unique cluster merger, 
that enable a direct detection of dark matter, independent of assumptions regarding the nature of the 
gravitational force law. Due to the collision of two clusters, the dissipationless stellar component and 
the fluid-like X-ray emitting plasma are spatially segregated. By using both wide-field ground based 
images and HST/ACS images of the cluster cores, we create gravitational lensing maps which show 
that the gravitational potential does not trace the plasma distribution, the dominant baryonic mass 
component, but rather approximately traces the distribution of galaxies. An 8a significance spatial 
offset of the center of the total mass from the center of the baryonic mass peaks cannot be explained 
with an alteration of the gravitational force law, and thus proves that the majority of the matter in 
the system is unseen. 

Subject headings: Gravitational lensing - Galaxies: clusters: individual: 1E0657-558 - dark matter 



1. INTRODUCTION 

We have known since 1937 that the gravitational po- 
tentials of galaxy clusters are too deep to be caused by 
the detected bar yonic mass an d a Newtonian gravita- 
tional force law l|Zwickvll937f) . Proposed solutions either 
invo ke dominan t quantities of non-luminous "dark mat- 
ter" l)OortF l9321 or alterations to either the gravitational 
force law |Bekcnstcin 2004; Brownstejn ig_Moffat 20QJ) 
or the particles' dynamical response to it l|Milgromll983|) . 
Previous works aimed at distinguishing between the 
dar k matter and alternative gravity hypothes es in galax- 
ies (|Buote ct al. 2002; Hockstra ct al. 2004") or galaxy 
clusters (Gavazzi 2002; Pointccoutcau & Silk 2005) have 
used objects in which the visible baryonic and hypoth- 
esized dark matter are spatially coincident, as in most 
of the Universe. These works favor the dark matter hy- 
pothesis, but their conclusions were necessarily based on 
non-trivial assumptions such as symmetry, the location 
of the center of mass of the system, and/or hydrostatic 
equilibrium, which left room for counterarguments. The 
actual existence of dark matter can only be confirmed 
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either by a laboratory detection or, in an astronomical 
context, by the discovery of a system in which the ob- 
served baryons and the inferred dark matter are spatially 
segregated. An ongoing galaxy cluster merger is such a 
system. 

Given sufficient time, galaxies (whos e stellar compo- 
nent makes up ~ 1 — 2% of the mass ijKochanek et alJ 
12003(1 under the assumption of Ne wtonian gra.vity) , 
plasma (~ 5 - 15% of the mass (|Allen et al.l l2002t 
iVikhlinin et alJl2006(l '). and any dark matter in a typical 
cluster acquire similar, centrally-symmetric spatial dis- 
tributions tracing the common gravitational potential. 
However, during a merger of two clusters, galaxies be- 
have as collisionless particles, while the fluid-like X-ray 
emitting intracluster plasma experiences ram pressure. 
Therefore, in the course of a cluster collision, galaxies 
spatially decouple from the plasma. We cl early see this 
effect in the unique cluster 1E0657-558 ((Tucker et alJ 
[iS9i). 

The cluster has two primary galaxy concentrations sep- 
arated by 0.72 Mpc on the sky, a less massive {T ^ 6 keV) 
western subcluster a nd a more massive (T ~ 14 keV) 
eastern main cluster llMarkevitch et aril2002() . Both con- 
centrations have associated X-ray emitting plasma off- 
set from the galaxies toward the center of the system. 
The X-ray image also shows a prominent bow shock on 
the western side of the western plasma cloud, indicat- 
ing that the subcluster is currently moving away from 
the main cluster at ~4700 km/s. As the line-of-sight ve- 
locity differe nce between t he components is only ~ 600 
km/s ( Ba rrena et al.ll2002() . the merger must be occur- 
ring nearly in the plane of the sky and the cores passed 
through each other ^ 100 Myr ago. 

Two galaxy concentrations that correspond to the 
main cluster and the smaller subcluster have moved 
ahead of their respective plasma clouds that have been 
slowed by ram pressure. This phenomenon provides an 
excellent setup for our simple test. In the absence of 
dark matter, the gravitational potential will trace the 
dominant visible matter component, which is the X-ray 
plasma. If, on the other hand, the mass is indeed domi- 
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Fig. 1. — Shown above in the top panel is a color image from the Magellan images of the merging cluster 1E0657— 558, with the white 
bar indicating 200 kpc at the distance of the cluster. In the bottom panel is a 500 ks Chandra image of the cluster. Shown in green contours 
in both panels are the weak lensing k reconstruction with the outer contour level at k = 0.16 and increasing in steps of 0.07. The white 
contours show the errors on the positions of the k peaks and correspond to 68.3%, 95.5%, and 99.7% confidence levels. The blue +s show 
the location of the centers used to measure the masses of the plasma clouds in Table 1^ 



nated by coUisionless dark matter, the potential will trace 
the distribution of that component, which is expected 
to be spatially coincident with the coUisionless galax- 
ies. Thus, by deriving a map of the gravitational po- 
tential, one can discriminate between these possibilities. 
We publishe d an initial attemp t at this using an archival 
VLT image l|Clowe et al.ll2004) : here we add three addi- 
tional optical image sets which allows us to increase the 
significance of the weak lensing results by more than a 
factor of 3. 

In this paper, we measure distances at the rcdshift of 
the cluster, z = 0.296, by assuming an fim = 0.3, A = 
0.7, Hq = 70km/s/Mpc cosmology which results in 4.413 
kpc/" plate-scale. None of the results of this paper are 
dependent on this assumption; changing the assumed 
cosmology will result in a change of the distances and 
absolute masses measured, but the relative masses of 
the various structures in each measurement remain un- 
changed. 

2. METHODOLOGY AND DATA 

We construct a map of the gravit ational poten- 
tial using wea k gravitational lensing ijMellierl Il999t 
IBartelmann &: Schnciderl l2001f) . which measures the dis- 
tortions of images of background galaxies caused by the 
gravitational deflection of light by the cluster's mass. 
This deflection stretches the image of the galaxy pref- 
erentially in the direction perpendicular to that of the 
cluster's center of mass. The imparted ellipticity is typi- 
cally comparable to or smaller than that intrinsic to the 
galaxy, and thus the distortion is only measurable statis- 
tically with large numbers of background galaxies. To do 
this measurement, we detect faint galaxies on deep op- 
tical images and calculate an ellipticity from the second 
moment of their surface brightness distribution, correct- 
ing the ellipticity for smearing by the point spread func- 
tion (corrections for both anisotropics and smearing are 
obtained using an im plementation of the KSB t echnique 
l|Kaiser et al.l Il995j) discussed in IClowe et al.1 (|Ip06)). 
The corrected ellipticities are a direct, but noisy, mea- 
surement of the reduced shear 3 = 7/(1 — k). The shear 
7 is the amount of anisotropic stretching of the galaxy 
image. The convergence k is the shape-independent in- 
crease in the size of the galaxy image. In Newtonian 



gravity, k is equal to the surface mass density of the lens 
divided by a scaling constant. In non-standard gravity 
models, k is no longer linearly related to the surface den- 
sity but is instead a non-local function that scales as the 
mass raised to a power less than one for a planar lens, 
reaching the limit of one half for constant acceleration 
ijMortlock fc Turneill200ll IZhao et alJl2006|) . While one 
can no longer directly obtain a map of the surface mass 
density using the distribution of k in non-standard grav- 
ity models, the locations of the k peaks, after adjusting 
for the extended wings, correspond to the locations of 
the surface mass density peaks. 

Our goal is thus to obtain a map of k. One can combine 
derivatives of g to obtain ( Schneider 1995; Kaiser 1995^) 

Vln(l-K) = - 1 52 \(9i,i+92,2\ 

1 - 5? - 5i V 52 1 - .91 A 52,1 - 51.2 ; 
which is integrated over the data fleld and converted into 
a two-dimensional map of k. The observationally un- 
constrained constant of integration, typically referred to 
as the "mass-sheet degeneracy," is effectively the true 
mean of ln(l — k) at the edge of the reconstruction. This 
method does, however, systematically underestimate k 
in the cores of massive clusters. This results in a slight 
increase to the centroiding errors of the peaks, and our 
measurements of k in the peaks of the components are 
only lower bounds. 

For 1E0657— 558, we have accumulated an exception- 
ally rich optical dataset, which we will use here to mea- 
sure g. It consists of the four sets of optical images shown 
i n_ Ta ble Q and the VLT image set used in Clowe et aO 
(2Q04.); the additional images signiflcantly increase the 
maximum resolution obtainable in the k reconstructions 
due to the increased number of background galaxies, 
particularly in the area covered by the ACS images, 
with which we measure the reduced shear. We reduce 
each image set independently and create galaxy cata- 
logs with 3 passband photometry. The one exception 
is the single passband HST pointing of main cluster, 
for which we measure colors from the Magellan images. 
Because it is not feasible to measure redshifts for all 
galaxies in the field, we select likely background galax- 
ies using magnitude and color cuts (m8i4 > 22 and not 
in the rhombus defined by 0.5 < mgoe — "^8i4 < 1-5, 
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TABLE 1 

Optical Imaging Sets 
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Note. — Limiting magnitudes for completion are given for galaxies and measured by where the number counts depart from 
a power law. All image sets had objects detected in the reddest passband available. 



'71435 - W606 > 1-5 * (mgoe - msii) - 0.25, and 771435 - 
"^■606 < l-6*(7n606 — "^8i4) + 0.4 for the ACS images; sim- 
ilar for the other image sets) that were cahbrated with 
photometric redshifts from the HDF-S ijFontana et al.l 
11999) . Each galaxy has a statistical weight based on its 
signif icance of detection in the image set owe et all 
12001 . and the weights are normalized among catalogs 
by comparing the rms reduced shear measured in a re- 
gion away from the cores of the cluster common to all five 
data sets. To combine the catalogs, we adopt a weighted 
average of the reduced shear measurements and appro- 
priately increase the statistical weight of galaxies that 
occur in more than one catalog. 

3. ANALYSIS 

We use the combined catalog to create a two- 
dimensional K reconstruction, the central portion of 
which is shown in Fig. ^ Two major peaks are clearly 
visible in the reconstruction, one centered 7'.'1 east and 
6'.'5 north of the subcluster's brightest cluster galaxy 
(BCG) and detected at 8 cr significance (as compared to 
3(7 in l)Clowe et al.l l200'3) ). and one centered 2'.'5 east 
and 11"5 south of the northern BCG in the main clus- 
ter (21'.'2 west and 17'.'7 north of the southern BCG) 
detected at 12a. We estimate centroid uncertainties by 
repeating bootstrap samplings of the background galaxy 
catalog, performing a k reconstruction with the resam- 
pled catalogs, and measuring the centroid of each peak. 
Both peaks are offset from their respective BCG by ~ 2a, 
but are within Icr of the luminosity centroid of the re- 
spective component's galaxies (both BCGs are slightly 
offset from the center of galaxy concentrations). Both 
peaks are also offset at ^ %a from the center of mass 
of their respective plasma clouds. They are skewed 
toward the plasma clouds, which is expected because 
the plasm a contributes about 1/lOth of the tota l clus- 
ter mass ()Allen et al.ll2002t iVikhlinin et alJl200fiD (and 
a higher fraction in non-standard gravity models with- 
out dark matter). The skew in each k peak toward the 
X-ray plasma is significant even after correcting for the 
overlapping wings of the other peak, and the degree of 
skewness is consistent with the X-ray plasma contribut- 
ing 14%^i[^^ of the observed k in the main cluster and 

10%^']pS in the subcluster (see Table|2|. Because of the 



large size of the reconstruction (34' or 9 Mpc on a side) , 
the change in k due to the mass-sheet degeneracy should 
be less than 1% and any systematic effects on the cen- 
troid and skewness of the peaks are much smaller than 
the measured error bars. 

The projected cluster galaxy stellar mass and plasma 
mass within 100 kpc apertures centered on the BCGs 
and X-ray plasma peaks are shown in Table El This 
aperture size was chosen as smaller apertures had sig- 
nificantly higher kappa measurement errors and larger 
apertures resulted in significant overlap of the apertures. 
Plasma masses were computed from a multicomponent 3- 
dimensional cluster model fit to the Chandra X-ray image 
(details of this fit will be given elsewhere). The emission 
in the Chandra energy band (mostly optically-thin ther- 
mal bremsstrahlung) is proportional to the square of the 
plasma density, with a small correction for the plasma 
temperature (also measured from the X-ray spectra), 
which gives the plasma mass. Because of the simplic- 
ity of this cluster's geometry, especially at the location 
of the subcluster, this mass estimate is quite robust (to 
a 10% accuracy). 

Stellar masses are calculated from the /-band lumi- 
nosity of all galaxies equal in brightness or fainter than 
the component BC G. The luminosities wer e converted 
into mass assuming ijKauffmann et alJl200^ M/Li = 2. 
The assumed mass-to-light ratio is highly uncertain (can 
vary between 0.5 and 3) and depends on the history of 
recent star formation of the galaxies in the apertures; 
however even in the case of an extreme deviation, the X- 
ray plasma is still the dominant baryonic component in 
all of the apertures. The quoted errors are only the errors 
on measuring the luminosity and do not include the un- 
certainty in the assumed mass-to-light ratio. Because we 
did not apply a color selection to the galaxies, these mea- 
surements are an upper limit on the stellar mass as they 
include contributions from galaxies not affiliated with the 
cluster. 

The mean k at each BCG was calculated by fitting a 
two peak model, each peak circularly symmetric, to the 
reconstruction and subtracting the contribution of the 
other peak at that distance. The mean k for each plasma 
cloud is the excess k after subtracting off the values for 
both peaks. 
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TABLE 2 

Component Masses 



Component 



RA (J2000) 



Dec (J2000) 



A/x(1O12M0) 



M*(1O^^M0) 



Main cluster BCG 
Main cluster plasma 
Subcluster BCG 
Subcluster plasma 



06 : 58 : 35.3 
06 : 58 : 30.2 
06 : 58 : 16.0 
06 : 58 : 21.2 



-55 : 56 ; 56.3 
-55 : 56 : 35.9 
-55 : 56 ; 35.1 
-55 : 56 ; 30.0 



5.5 ±0.6 

6.6 ±0.7 

2.7 ±0.3 

5.8 ±0.6 



0.54 ±0.08 
0.23 ± 0.02 
0.58 ±0.09 
0.12 ±0.01 



0.36 ±0.06 
0.05 ±0.06 
0.20 ±0.05 
0.02 ±0.06 



Note. — All values are calculated by averaging over an aperture of 100 kpc radius around the given position (marked with 
blue +s for the centers of the plasma clouds in Fig0. R measurements for the plasma clouds are the residual left over after 
subtraction of circularly symmetric profiles centered on the BCGs. 



The total of the two visible mass components of the 
subcluster is greater by a factor of 2 at the plasma peak 
than at the BCG; however, the center of the lensing mass 
is located near the BCG. The difference of the baryonic 
mass between these two positions would be even greater 
if we excluded a contribution of the non-peaked plasma 
component between the shock front and the subcluster. 
For the main cluster, we see the same effect, although 
the baryonic mass difference is smaller. Note that both 
the plasma mass and the stellar mass are determined 
directly from the X-ray and optical images, respectively, 
independently of any gravity or dark matter models. 

4. DISCUSSION 

A key limitation of the gravitational lensing method- 
ology is that it only produces a two-dimensional map of 
K, and hence raises the possibility that structures seen in 
the map are caused by physically unrelated masses along 
the line-of-sight. Because the background galaxies reside 
at a mean z ~ 1, structures capable of providing a signifi- 
cance amount of k must lie at z < 0.8. By comparing the 
measured shear for galaxies divi ded into crude redshif t 
bins using photometric redshifts (|Wittman et a,].ll2nnl . 
we further limit the redshift of the lensing objects to 
0.18 < z < 0.39. This range is consistent with the clus- 
ter redshift, but corresponds to a large volume in which 
a structure unassociated with the cluster could exist and 
be projected onto the lensing map. However, the num- 
ber density of structures with these lensing strengths in 
blank field surveys is ~ 10"'^ arcmin"^ for the subcluster 
fWittman ct al. 2003), and an order of magnitude less 
for the main cluster, resulting in a 10~^ probability of 
having two structures within a square arcminute of the 
observed cluster cores. Further, all such lenses observed 
in these cosmic shear surveys are clusters with enough 
plasma and galaxies to be easily observable. There is 
no evidence, however, in our deep imaging for additional 
cluster sized concentrations of galaxies or of plasma hot- 
ter than r 0.5 keV (the lower bound of the Chandra 
energy band) near the observed lensing structures. 

Another alternate explanation of the lensing signal is 
related to the fact t hat clusters form a t the intersections 
of matter filaments ijBond et alJIlQQ^ . In principle, one 
could imagine two line-of-sight filaments of intergalactic 
gas (too cool to be visible with Chandra and too dif- 
fuse to have cooled into stars) extending from the clus- 
ter at the locations of the weak lensing peaks. To ex- 
plain the measured surface mass density, such filaments 
would have to be several Megaparsecs long, very nar- 



row, and oriented exactly along the line of sight. The 
probability of such an orientation for two such filaments 
in the field is ^ 10~^. Further, because the two clus- 
ter components are moving at a relative transverse ve- 
locity of 4700 km/s compared to the typical peculiar 
velocities in the CMB frame of a few hundred km/s, 
the filaments could coincide so exactly with each of the 
BCGs only by chance. This is an additional factor of 
~ 10~^ reduction in probability. While such projections 
become more important in non-standard gravity mod- 
els because in such models the thin len s approximation 
breaks down l|Mortlock fc Turner! I2001D and structures 
with a given surface density produce a greater amount 
of lensing the more they are extended along the line-of- 
sight, two such projections would still have a <C 10~* 
probability. Finally, we r nention th at two other merg- 
ing clusters, MS1054-03 ijJee et al.l l200^ and A520 (in 
preparation), exhibit similar offsets between the peaks of 
the lensing and baryonic mass, although based on lens- 
ing reconstructions with lower spatial resolution and less 
clear-cut cluster geometry. 

A final possibility is that some alternative gravity mod- 
els may be able to suppress the lensing potential of the 
centra l peak in a multiple-peak system, as in Angus ct aQ 
( 2006) . That work used a model of a gas disk located be- 
tween two symmetric mass concentrations representing 
the galaxy subclus ters. In their k map, derived in the 
TeVeS framework l)Bekensteinll2004j) . the relative signal 
from this disk may be suppressed, but would still be eas- 
ily visible with the noise levels of our reconstruction. Our 
K map, however, has no evidence of any mass concentra- 
tion between the two galaxy subclusters other than the 
small perturbations consistent with the gas mass contri- 
bution in Newtonian gravity. Furthermore, such a sup- 
pression has also only been shown to work for symmetric 
systems which have the central peak directly between 
the two outer peaks. In 1E0657— 558, however, the X- 
ray plasma, which would provide the central peak, lies 
north of the line connecting the two k peaks. Further, the 
absolute k levels of the peaks observed in 1E0657— 558 
are in good agreement with those in systems with sim- 
ilar velocity dispersions and X-ray temperatures (e.g. 
iClowe fc Schneider! I2002D which have the gas and the 
galaxies coincident. The K-to-light ratios are also con- 
sistent with those in normal clusters with coincident gas 
and galaxies. Therefore one would need to not only sup- 
press the inner peak in the k map relative to the two 
outer peaks in this system, but also enhance the strength 
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of the outer peaks to make up for the missing plasma 
mass. 

Any non-standard gravitational force that scales with 
baryonic mass will fail to reproduce these observations. 
The lensing peaks require unseen matter concentrations 
that are more massive than and offset from the plasma. 
While the existence of dark matter removes the pri- 
mary motivation for alternative gravity models, it does 
not preclude non-standard gravity. The scaling relation 
between n and surface mass density, however, has im- 
portant consequences for models that mix dark matter 
with non-Newtonian gravity: to achieve the ~ 7 : 1 
ratio in n between the dark matter -I- galaxy compo- 
nent and the plasma component (Table \^ , the true ra- 
tio of mass would be even higher (as high as 49:1 for a 
constant acceleration model; although MOND l)Milgroml 
[TqSS) would not reach this ratio as the dark matter 
density would become high enough to shift the accel- 
eration into the quasi-Newtonian regime), making the 
need for dark matter even more acute. Such high con- 
centrations of dark matter, however, are extremely un- 
likely based on the meas ured X-ray plasma temperatures 
l)Markevitch et ani2002D and cluster galaxy velocity dis- 



persions l|Barrena et al.ll2002|) . 

The spatial separation of the dominant baryonic com- 
ponent in a galaxy cluster from the hypothesized dark 
matter produced during a cluster merger has enabled 
us to directly compare the dark matter hypothesis to 
one with only visible matter but a modified law of grav- 
ity. The observed displacement between the bulk of the 
baryons and the gravitational potential proves the pres- 
ence of dark matter for the most general assumptions 
regarding the behavior of gravity. 
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